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Abstract 

High density magnetized detectors are well suited to exploit the outstanding 
purity and intensities of novel neutrino sources like Neutrino Factories and Beta 
Beams. They can also provide independent measurements of leptonic mixing pa- 
rameters through the observation of atmospheric muon-neutrinos. In this paper, 
we discuss the combination of these observables from a multi-kton iron detector 
and a high energy Beta Beam; in particular, we demonstrate that even with 
moderate detector granularities the neutrino mass hierarchy can be determined 
for #13 values greater than 4°. 



PACS: 14.60.Pq, 14.60.Lm 



1 Introduction 



In the standard interpretation of the experimental evidence for neutrino oscillation [1, 
12], the squared- mass differences among the v mass eigenstates (m 1; m 2 and m?) are 
rather hierarchical: Am| x = m\ — m\ ~ 1/30 x |Am3 2 | ~ jAm^l- The small size 
of Amjj is implied by solar and long-baseline reactor neutrino data. In this frame- 
work, particularly enlightening are the oscillations at the "atmospheric scale", i.e. 
when \Amf 2 \L/4E ~ tt/2, L and E being the neutrino energy and path- length, respec- 
tively. A muon-neutrino oscillating at the atmospheric scale undergoes mainly — + v T 
transitions. This dominant mode is implied by long-baseline accelerator experiments, 
atmospherics and short-baseline reactor data; moreover, it is currently under test in a 
direct manner at CNGS [3]. We have no evidence of subdominant — > v e transitions 
at such scale, implying a small mixing angle between the first and third mass eigenstate 
(#13 < 10° at 90% C.L. [1]). If #13 is non-zero, the subdominant — > v e amplitude 
and its T or CP conjugates encode a wealth of information. In particular, it allows the 
determination of #13 and the Dirac complex phase of the leptonic mixing matrix. The 
— > v e transition probability is also perturbed by matter effects if the path of the 
neutrinos through the Earth is sufficiently large. The perturbation depends on the sign 
of Am|i and therefore it allows the determination of the hierarchy among the neutrino 
masses. These considerations explain the enormous interest toward novel neutrino 
sources operating at the atmospheric scale [S] like, for instance, Neutrino Factories [5] 
or Beta Beams [7J. Combination of these facilities with more traditional measurements 
of atmospheric neutrinos has also been considered [8], [9], [TO] . Such combination is par- 
ticularly natural when the detector at the far location is dense and capable of charge 
reconstruction: this is the case of magnetized iron calorimeters, which can perform 
detailed measurements of atmospheric fluxes and are recognized as ideally suited to 
exploit a possible Neutrino Factory or high energy Beta Beam [HI [121 E31 EH E5] ■ How- 
ever, it is generally believed that the contribution of atmospheric neutrinos is marginal 
for any realistic configuration, unless the value of the #13 angle turns out to be very 
close to current bounds [16J. In this paper, we carry out a more detailed analysis of 
the atmospheric data that can be collected by a 40-kton magnetized iron calorimeter 
(SecfS]). This analysis is tuned to identify the occurrence of resonant transitions in 
the earth and it is combined with the data that can be collected by the same detec- 
tor exposed to a high energy Beta Beam [T2] at a baseline of ~700 km (SecHJ). The 
sensitivity of the atmospheric data to the neutrino mass hierarchy does not depend on 
the value of 5 ( "CP-blind" ) but it is highly deteriorated by the lack of knowledge of 
#13; on the other hand, the Beta Beam data exhibits a very strong dependence on 5 
but provide tight constraints on the size of the #13. Their combination (Sec. [5]) results 
in a significantly improved capability to determine the hierarchy of the neutrino mass 
eigenstates. 
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2 Atmospheric and long-baseline oscillations 



Current experimental data are unable to fix uniquely the hierarchy of neutrino masses pQ, 
i.e. whether the mi eigenstate is lighter than m 3 (mi < m 2 < m 3 : "normal hierarchy") 
or heavier (m 3 < mi < m 2 : "inverted hierarchy"). Clearly, this pattern is of great the- 
oretical relevance because it allows discrimination among models that explain neutrino 
masses [T7j. In vacuum, observables sensitive to the mass pattern exist but their ex- 
ploitation is extremely challenging from the experimental point of view [18] . However, 
propagation in matter can enhance the perturbations to the transition probabilities due 
to the sign of Am^ and, for sizable values of 6*i 3 such perturbations become observable 
even with present technologies. In particular, the v e — > transition amplitude (or its 
T and CP-conjugate) encodes an explicit dependence on the sign of Am 31 [19, 20|: 
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In this formula A = Am^L / (AE) and the terms contributing to the Jarlskog invariant 
are split into the small parameter sin2#i 3 , the 0(1) term £ = cos(9i 3 sin2#i2 sin 2023 
and the CP term sin 6; A = 2V2G F n e E/Am 2 31 with Gp the Fermi coupling constant 
and n e the electron density in matter. Note that the sign of A depends on the sign of 
Arn^i which is positive (negative) for normal (inverted) hierarchy of neutrino masses. 

For a monochromatic beam, a v e — > z/ M run combined with its CP conjugate v e — > 
cannot determine uniquely sign(Am| 1 ) for all values of 5. In particular, if the hierarchy 
is normal (inverted) and 5 < (5 > 0), it is always possible to find a solution that 
reproduces the correct v e — > and v e — > rate assuming the wrong hypothesis on 
sign(Am 31 ). For beams of finite width, spectral information help lifting this ambiguity 
but, in general, they require very large statistics and excellent detector resolution. 
Figs. [T] and [2] show the allowed #i 3 , 5 region from a neutrino and anti- neutrino run of 
a Super-SPS [21] based Beta Beam assuming both the wrong (red) and right (green) 
hypothesis^- The true parameters are #i 3 = 4°, sign(Am| 1 ) = +1 and S = 90° (FigfT]) or 
5 = —90° (Fig E]). The existence of overlapping regions in Fig ]2] testifies the consistency 
of the wrong hypothesis with the accelerator data. In general, we expect a highly 

1 Details on this experimental setup are given in [12] and briefly recalled in Sec. [3l 
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deteriorated sensitivity to mass hierarchy for nearly 50% of the possible true values of 
5. This is clearly visible in Fig. 8 of [12] and it is reproduced in Fig 16] (blue linej§. 
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Figure 1: Allowed #13, 5 regions from a neutrino and anti-neutrino run of a Super-SPS 
based Beta Beam assuming the wrong (inverted - red) and right (normal - green) mass 
hierarchy hypothesis. The true parameters are # 13 = 4°, sign(Am| 1 ) = +1 and 8 = 90°. 
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Figure 2: As in Fig. [T] with 8 
wrong hypothesis (red area). 



-90°. Overlapping regions are present also for the 



Since the dependence on 8 and, in general, the three-family interference effects are 
the origin of this ambiguity, a substantial improvement can be achieved exploiting CP- 



2 Actually, most of the next generation long-baseline experiments proposed so far suffer from this 
limitation. For an example based on Superbeams we refer to [22] . 
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blind observables. In u e — ► transitions, CP-blindness can be achieved introducing 
a phase advance (increase of baseline) such that the phase acquired by neutrinos due 
to interaction with matter equals 2tt [23]. It can be shown [24] that this condition 
is equivalent to choosing L = 2ix j \f2Gpn e , i.e. finding a L such that the last three 
terms of EqJT] cancels. Such baseline is dubbed "magic" in literature [23] and it is 
proportional to the energy and width of the MSW resonance in a medium of constant 
density: 

E R = ±Am 2 31 L 

magic 

COS 2^13/ 4:71 (1) 

T R = \Amn\Lmagic sm29 13 /2n (2) 

the + (— ) sign referring to neutrinos (antineutrinos). A similar approach can be pur- 
sued using atmospheric neutrinos. In this case, however, additional difficulties are 
present. First of all, the initial flux is a mixture of v^, v e and their antiparticles. More- 
over, the approximation of constant density is non tenable and core-mantle interference 
effects are possible. Finally, the detector resolution must be sufficient to identify the 
region where the MSW resonance occurs. In turn, the latter condition is a require- 
ment on the size of #13 through Eqj2j Eq. [1] implies that the resonance condition can 
occur only for neutrinos if the hierarchy is normal and only for anti-neutrinos if the 
hierarchy is inverted. Hence, the task of fixing the hierarchy is reduced to the issue 
of determining the occurrence of this resonance or at least the spectral perturbation 
induced by it in the neutrino (antineutrino) sample. Such perturbation is depicted in 
FigJSJfor sin 2 2^13 = 0.1 (normal hierarchy) as a function of L/E and the Nadir angle 
9 of the neutrinos. It has been computed solving numerically [25] the Schroedinger 
equation for the propagation of neutrinos in matter and assuming the Earth density 
profile of Ref. [26]. For all currently allowed values of Am^/I Amj^l, three family in- 
terference effects are negligible in the multi-GeV neutrino sample [IE] (E > 1.5 GeV) 
and, therefore, the observables are completely CP-blind. 



3 Experimental setup 

The experimental setup considered in this paper closely follows the configuration stud- 
ied in Ref. [12]. It exploits a possible novel source operating in v e —>■ mode, i.e. a 
Beta Beam leveraging an upgraded injection complex for the LHC. The Beta Beams 
are pure sources of electron-neutrinos obtained producing, accelerating and stacking 
beta-unstable isotopes. The most advanced machine study regarding this novel ap- 
proach is currently in progress within the Eurisol Design Study group [27]; it is based 
on existing CERN accelerating machines (PS and SPS) and uses 6 He and 18 Ne as v e 
and v e emitters, respectively. 

Compared with such baseline configuration, the setup of [12] . which was originally 
proposed in Ref. [11] for a massive water Cherenkov detector, considers a higher energy 
proton injector to the LHC (the "Super-SPS" instead of the existing SPS) employed 
to increase the 7 of the ions; in this case, the neutrinos have energies exceeding 1 GeV. 
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Figure 3: Probability of observing atmospheric muon neutrinos (left plot - antineutrinos 
are shown in the right plot) from — > and v e — > z/ M transitions as a function of L/E 
and the Nadir angle 9 (sin 2 2^13 = 0.1, Am| x = 2.5 • 10~ 3 eV 2 ). In the upper label, r 
is the ratio between the initial and v e fluxes. 



The corresponding spectra are recalled in Fig. HI They are computed assuming the far 
detector to be located 730 km from CERN (CERN-to-Gran Sasso distance) and the ions 
accelerated up the same 7 (7 = 350) both for 6 He and for 18 Ne ("7 = 350, 350 option"[§) 
or up to the maximum rigidity allowed by the Super-SPS (7 = 580 for 18 Ne, 7 = 350 for 
6 He) . At these energies, final state muons produced by charged-current interactions 
have a range significantly larger than the pion interaction length in iron. Hence, muon 
identification against the bulk of v e CC and NC interactions becomes possible even 
for high density iron detectors. Differently from the case of the Neutrino Factories, 
the far detector does not require charge identification capability: the initial state v e 
beam is free from z/ M contaminations and n + production from charm or tau decays is 
kinematically suppressed. In this context, the only advantage of magnetization comes 
from the additional suppression of the NC background with a pion misidentified as a 
muon. 

In Ref. [12], the detector was based on glass Resistive Plate Chambers interleaved 
with 4 cm thick vertical iron slabs. The fiducial mass is 40 kton i.e. the detector 
can be accommodated in one of the existing LNGS experimental halls. The RPC 
are housed in a 2 cm gap and the signal is readout on external pick-up electrodes 
segmented in 2 x 2 cm 2 pads. This detector can be divided into large modules (about 
15 m length) and magnetized by copper coils running through the slabs. This design, 
which resembles the toroidal configuration of MINOS, was studied in detail by the 

36 He and 18 Ne having nearly the same Q-value, equal 7 corresponds to equal average neutrino 
energy. 
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Figure 4: Neutrino spectra at the far location. 



MONOLITH Collaboration in 1999 [28J and it has been recently revived by the Indian 
proposal INO j29j[30]. Hence, the far detector of [32] represents a unique atmospheric 
neutrino detector. In particular, MONOLITH studied have demonstrated that the 
vertical orientation of the steel does not compromise the capability of reconstructing 
the oscillation pattern of atmospheric neutrinos. Moreover, the atmospheric analysis 
profits of the higher granularity needed for the accelerator-based appearance search. 

In the following, we consider a 40 kton magnetized detector running* at a Beta 
Beam for 10 years. The accelerator facility is assumed to provide 2.9 xlO 18 6 He and 
l.lxlO 18 18 Ne decays per year. As in Ref. [12], we stress that no solid estimate of 
fluxes are available for the Super-SPS option: we refer to [L2] for a systematic study 
of the sensitivity as a function of fluxes. 



4 Analysis of atmospheric data 

The simulation of atmospheric events at the far detector for the Beta Beam is similar 
to the one recently implemented by MINOS [31]: flux expectations are based on the 
Bartol 96 [32] model and interactions have been calculated with GRV94 [33] parton 
distributions including contributions from quasi-elastic and single pion production [25] . 
The detector response has been fully simulated with GEANT3 [34] assuming the mag- 
netic field maps resulting from the coil arrangement of MONOLITH. The average 
magnetic field along the slabs is ~1.3 T. 

4 The interest for a vertical slab configuration was originally motivated by the possible exploitation 
of MONOLITH at CNGS (see Ref. [25] Chap.2). 
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As discussed in Sec. [21 a substantial enhancement of the sensitivity to hierarchy can 
be obtained tuning the event selection to identify the region where resonant enhance- 
ment occurs. The resonance appears as a perturbation to the original L/E pattern for 
a given earth matter density. For cos# > 0.85 the neutrinos cross the earth core and 
resonant conversions can appear already at few GeV (En ~ 3 GeV). At smaller cos 6 
v's intersects the Earth mantle and crust and resonances appear at larger Er (~ 7 
and ~ 11 GeV respectively). Hence, a two dimensional analysis either in E versus 
L or, equivalently, in L/E versus cos6 is appropriate to extract information on the 
existence and location of Earth matter resonances. The choice of L/E a variable al- 
lows the implementation of the Monolith criteria for the observability of the sinusoidal 
pattern [28] to the search for resonant conversion. This approach has been developed 
in [25], where the expected resolution on L and E is computed on an event-by-event 
basis; here, events are retained only if the expected resolution is adequate (FWHM 
on L/E smaller than 50%) in the region of interest for matter effect perturbations. 
A preselection is applied requiring a minimum muon energy of 1.5 GeV; only internal 
events are considered: they can be either fully contained or with a single outgoing track 
ranging out of the detector with a visible path-length greater than 4 m. On average, 
these requirements result in a selection efficiency that is negligible below 3 GeV and 
nearly constant (~55%) above 5 GeV. At these energies, charge misidentification is 
well below 5% and therefore, rather pure /i + and pT samples can be collected. 

Fig. [5] (left) shows the L/E distribution of selected events corresponding to 400 
kton-year (10 y of data taking) for ISm\\ = 2.5- 10~ 3 eV 2 , sin 2 26 13 = 0.1 and 9 23 = 45°. 
In this figure, also down-going events are shown: they are used as an unoscillated refer- 
ence sample to reduce the systematics due to the uncertainty on the initial fluxes. The 
ratio upgoing/downgoing is shown in Fig. [5] (right); the fit assumes 6*13 = 0°: pertur- 
bations due to matter effects and finite #13 are visible only for neutrinos since normal 
hierarchy is assumed; moreover, distortion with respect to pure — > v T oscillations 
are present mainly beyond the first minimum of the survival probability. 

5 Sensitivity to mass hierarchy 

The determination of mass hierarchy from atmospheric neutrinos is plagued by the 
uncertainty on the mixing parameters, particularly on # 13 . Accelerator neutrino ex- 
periments (a Beta Beam in the present case) can provide very strong constraints on 
this angle, once marginalized over 5. In particular, the manifold of the allowed #13, S 
values is highly reduced by running both in u e — > and in its CP-conjugate mode 
(anti- neutrino run). Hence the combination of accelerator and atmospheric data can 
be rather effective 0, HO] . 
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Figure 5: L/E distribution of selected events corresponding to 400 kton-year (10 y of 
data taking) for Am^ = 2.5 ■ 10~ 3 eV 2 , sin 2 29 13 = 0.1 and 6» 23 = 45° (left upper: 
neutrinos, left lower: antineutrinos) . The shaded (unshaded) histogram corresponds to 
the sample of up-going (down-going) neutrinos. For illustration of the effect, the ratio 
upgoing/downgoing is also shown (right) together with a fit that assumes #i 3 = 0°. 

5.1 Likelihood analysis 

The analysis of accelerator data employed in the following does not differ from the 
one described in |12j : we refer to that reference for the details of event selection and 
detector performance. The sensitivity to the sign of Am| x is computed from a likelihood 
analysis that combines the allowed #13 range resulting from the Beta Beam run with 
the atmospheric likelihood: 

ln£ = lnf] 

i,g 

where the subscripts i and q are referred to bins in the (cos#, L/E) plane and of muon 
charge respectively; U iiq and A q \i iA are the observed and expected number of up-going 
neutrino events in the i, q-th bin and the last term in the likelihood accounts for the 
overall normalization uncertainty for neutrinos and anti-neutrinos separately: further 
details on the treatment of systematics and their relevance for the mass hierarchy 
determination are provided below (Sec. 15. 2p . Additional penalty terms of the form 
(P — Pbest) 2 /2cx 2 are added to ([3]) to account for the uncertainty on P = [Am^] and 
P = #23- Here, Pbest is the current best estimate of the parameter and a the estimated 
precision. The latter is the current precision for #23 and the expected one at the end 
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of T2K for lAmjjJ. For each "true" value of #13, 5 and the sign of Am^ (e.g. +1 
for normal hierarchy), we compute the CL. of the best fit for the wrong hypothesis 
(e.g. sign(Am| 1 ) = —1). Fig. [6] shows the region where the wrong hypothesis is 
excluded at 90, 95 and 99% CL. This figure demonstrates what stated qualitatively in 
Sec. [2J due to the insensitivity to the Dirac phase, the cancellation effect that equals 
the rates of /i + and \i~ at the Beta Beam does not take place with atmospherics. 
Hence, atmospheric data are relevant and, actually, dominate the sensitivity to mass 
hierarchy in most of the Beta Beam blind region. In the present scenario, hierarchy 
can be determined at 90% CL. down to #13 ~ 4°. At positive (negative) S and normal 
(inverted) hierarchy, the Beta Beam data dominates and sign(Amg 1 ) can be extracted 
down to 9 1 3 ~ 2°. Compared with stand-alone analysis of atmospheric neutrinos with 
high density detectors [25], where some sensitivity can be achieved only for 13 values 
very close to current CHOOZ limits, here the accelerator-based determination of this 
parameter pushes the sensitivity to the mass hierarchy down to sin 2 2# 13 ~ 0.02 even 
for very unfavorable values of the Dirac phase 5. This is illustrated in Fig. [7J where the 
likelihoods for the atmospheric sample only are shown for # 13 = 5°. The left (right) plots 
assume normal (inverted) hierarchy as the true hierarchy. In general, Fig. [7J shows that 
hierarchy separation cannot be achieved if #13 is completely unconstrained: on the other 
hand, even for the most unfavorable values of S, the Beta Beam provides stringent limits 
to #i3 that makes the atmospheric sample a sensitive tool for hierarchy determination. 
If the true hierarchy is inverted, distortions will appear in the antineutrino sample, 
i.e. in the atmospheric sample with lower statistics due to the difference between 
and Up cross sections. However, the discrimination power depends on the difference 
between the higher statistics sample and the lower statistics one and, as visible 
from Fig. [7J the power of the test is practically the same both assuming normal and 
inverted "true" hierarchy. 

The use of the Monolith criteria for event selection improves substantially the sen- 
sitivity with respect to analyses where flat resolutions have been assumed [16J. In spite 
of their limited mass, high density atmospheric detectors compete with megaton-size 
water Cherenkov's thanks to the capability of measuring the sign of the multi-GeV 
sample. On the other hand, due to larger statistics, accelerator data are sensitive to 
perturbations induced by the sign of Am 2 up to sin 2 26* 13 ~ 0.03 — 0.02 (see Fig. 16 of 
Ref.|10j). Finally, it is worth noting that significantly better results (sensitivities to the 
hierarchy down to sin 2 26>i 3 < 0.01) could be achieved by dedicated detectors located 
at the magic baseline without the use of atmospheric sample [14|. These setups, how- 
ever, requires either large boosters for the beta-beam ions (LHC) or novel ions with 
larger Q- values and increased flux to compensate for the losses at the far location [T5] . 
Moreover, no information on CP violation in the leptonic sector would be accessible 
and a second detector at a closer distance from the source would be needed. 
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Figure 6: Region of the true value of the parameters #13 and S where the correct 
neutrino hierarchy (assumed to be normal in this figure) can be distinguished from the 
wrong one at 90, 95 and 99% CL. Data from the Beta Beam only are also shown. 
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Figure 7: Likelihoods of the atmospheric sample (A% 2 = — 2Log(L/ L max )) as a function 
of #13 assuming normal hierarchy (left plot) or inverted hierarchy (right plot) as true 
hierarchy and #13 = 5°. The continuous green line is the likelihood for Am 2 > (right 
hypotesis in the left plot, wrong hypothesis in the right plot) and the dashed red line 
is for Am 2 < 0. 
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5.2 Systematic uncertainties for the atmospheric sample 

Systematic uncertainties on the rate of multi-GeV and sub-GeV samples of atmospheric 
neutrinos have been studied in depth in literature due to their relevance for the de- 
termination of the leading oscillation parameters Am 13 and sin 2 #23 . Moreover, they 
affect the pattern of muon disappearance and the perturbation of the v e sample and, 
therefore, they can be limiting factors in the determination of the mass hierarchy and 
#13. In our case, the knowledge of #13 is completely dominated by accelerator data 
while the determination of the mass hierarchy is mainly related to the location of the 
MSW resonance in the neutrino or antineutrino sample. Hence, several systematic 
contributions that currently limit the knowledge of the leading parameters at atmo- 
spheric detectors (overall flux and cross section normalization, v^j u e ratio, zenith and 
energy dependence of the differential fluxes) play here a minor role. In order to probe 
quantitatively such statement, in the present analysis systematics uncertainties have 
been treated following the approach developed by the Monolith collaboration and based 
on up-down differential normalization [28, 35]. For each bin in the L/E distribution 
of the up-going neutrinos, the unoscillated sample is taken from the corresponding 
down-going. Here "up-going" means that neutrinos reach the detector coming from 
below the horizon (cos 8 > 0); the "corresponding down-going" are neutrinos whose 
Nadir angle 9' is such that 9' = 71 — 9. This subsample provides a nearly theory-free 
normalization since the uncertainties on the flux and cross sections affecting the pop- 
ulation of such bin cancel out because of the spherical symmetry and the isotropy of 
the primary cosmic ray flux. Further details and interesting subtleties can be found 
in the Monolith proposal [28]. For instance, geomagnetic effects can introduce an 
asymmetry in the up/down sample that could be interpreted as a perturbation of the 
L/E pattern. However, for the multi-GeV sample selected in this analysis, with a 
minimum neutrino energy of 3 GeV (see Sec. HJ), this effect is negligible [36]. Yet, an 
overall penalty term in the flux normalization has been introduced (see Eq. |3J , to cope 
with the Up/ v e ratio uncertainty, as the proposed detector cannot directly measure the 
v e flux. This approach is somewhat conservative: in most of the cases, the present 
theoretical knowledge of the atmospheric fluxes exceed the statistical precision of the 
down-going samples for bin widths comparable to the detector resolution. However, 
for the present case where systematics play a less important role than for the leading 
parameter analysis, this technique simplifies remarkably the treatment of systematics 
without compromising the sensitivity to the mass hierarchy. To ease comparison with 
other analyses (See e.g. Table 1 of [16]) we note that the statistical uncertainty of 
the down-going sample in the region where the oscillation dip at #13 is located is of 
about 9%, while typical differential uncertainties in the rate of multi-GeV are of the 
order of 5%. Moreover, the down-going normalization sample is measured for neutrinos 
and antineutrinos separately, hence assuming that the uncertainty on the v^j ratio 
could also have a dependence in energy and in direction, while in most of the analyses 
a flat uncertainty of about 5% is considered. As a consistency check, we artificially 
decreased the statistical uncertainty of the down-going sample up to 3% noting that 
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the sensitivity to the sign of Am^ remains practically unchanged. This is in agree- 
ment to what already observed by some authors, i.e. that "the discrimination between 
normal and inverted hierarchy is based on a very characteristic signal consisting of 
pronounced structures in the E and cos 6 distributions, which cannot easily mimicked 
by the systematic effects" [16J. 

The most relevant systematic uncertainty of instrumental origin is related to the 
pollution of the neutrino sample with antineutrinos due to charge misassignment. In 
most of the analyses presented in literature, an overall uncertainty on the knowledge 
of the charge (typically of 5%) is assumed independently of the topology of event. 
The present analysis, instead, is based on a full GEANT3 simulation of the detector, 
followed by a reconstruction stage with track fitting. Neutrino flavour assignment is 
based on the reconstructed muon charge, thus encompassing a realistic energy and event 
shape dependency, as far as GEANT3 can predict the behaviour of magnetised iron 
calorimeter^]. The charge misidentification systematics is the most relevant detector- 
dependent systematics in the atmospheric analysis, while it is immaterial in the Beta 
Beam analysis, where ir/fi misidentification dominates. This is the reason why we 
assumed the systematics of the accelerator and atmospheric data fully uncorrelated. 



6 Conclusions 

Detailed atmospheric neutrino studies can be done in several detectors that have been 
built or proposed for long-baseline accelerator facilities. In particular, high density 
magnetized detectors can study multi-GeV atmospheric neutrinos and distinguish be- 
tween and Vp. These data can be combined in a highly non trivial manner to 
improve the sensitivity to the neutrino mass hierarchy. In this paper we considered 
in particular a high energy Beta Beam with a 40 kton iron calorimeter located at the 
CERN to Gran Sasso distance. The atmospheric data collected by this detector allows 
the determination of the neutrino hierarchy down to #13 ~ 4° in the region of 5 where 
the corresponding Beta Beam data cannot constrain the sign of Am^. Far from this 
blind region, the Beta Beam data dominates and sign(Amg 1 ) can be extracted down 
to #13 ~ 2° at 90% C.L. 
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